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The effects of vertical outflows on disk dynamos 
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Abstract. We consider the effect of vertical outflows on the mean-field dynamo in a thin disk. These outflows could 
be due to winds or magnetic buoyancy. We analyse both two-dimensional finite-difference numerical solutions of 
the axisymmetric dynamo equations and a free-decay mode expansion using the thin-disk approximation. Contrary 
to expectations, a vertical velocity can enhance dynamo action, provided the velocity is not too strong. In the 
nonlinear regime this can lead to super-exponential growth of the magnetic field. 
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1. Introduction 

Magnetic fields play an important role in accretion disks. 
Inside the disk, for example, they allow the magneto- 
rotational instability to develop ( Balbus fc Hawlcy 1998| ). 
This instability can lead to the development of turbu- 
lence which in turn transports angular momentum radi- 
ally outwards, thus allowing the matter to accrete. On 
the other hand, large scale magnetic fields are often con- 
sidered necessary for launching winds or jets from accre- 
tion disks. The origin of such large scale magnetic fields is 
still unclear: they may be advected from the surround- 
ing medium or be generated by a dynamo inside the 
disk. Advection appears unlikely as turbulence leads to 
enhanced viscosity and magnetic diffusivity, so that the 
two are of the same order (Prandtl number of order unity; 
Pouquct et al. 1976 ). Thus, turbulent magnetic diffusivity 
can compensate the dragging of the field by viscously in- 



duce d accretion flow (van Ballcgooijen 1989 : Lubow et al 



1994; Heyvaerts et al. 1996). Dynamo action is a plausi- 



ble mechanism for producing magnetic fields in accretion 
disks ( |Pudritz 198l| ; |Stepinski fc Levy 1988| ). However, 



dynamo magnetic fields can affect winds from accretion 
disks ( Campbell et al. 1998] ; |Campbcll 1999 ; Brandenburg 
et al. 2000). On the other hand, the wind can also af- 



fect the dynamo. In particular, an enhancement by wind 



was suggested by Brandenburg et al. (1993) in the con- 



text of galactic dynamos. In the same context, Elstncr 



et al. (1995) also consider the effects of shear in the wind. 

In this paper, we consider the generation of magnetic 
field by a dynamo acting in a thin accretion disk. We pro- 
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vide a study of the effects of vertical velocities. We show 
how vertical velocities can enhance the dynamo, allow- 
ing a larger growth rate and leading to super-exponential 
growth of the magnetic field in the nonlinear regime. 

The vertical velocities considered can have several ori- 
gins. They can, for example, be due to a wind emanating 
from the disk or to magnetic buoyancy inside the disk. 
Note that we are interested here in magnetic buoyancy 
as generating v ertical outflows. W e neglect the associated 
dynamo effect ( Moss et al. 1999 ). 

The paper is organized as follows. In Sect. 2, we de- 
scribe the two complementary approaches used and spec- 
ify our choice of parameters. In Sect. 3, we discuss the 
linear results without vertical outflows. Results with ver- 
tical outflows are discussed in Sect. 4 for the linear regime 
and in Sect. 5 for the nonlinear regime. Sect. 6 presents 
our conclusions. 



2. The models 

Simulations of dynamo-generated turbulence in accretion 
disks have shown that, if stratification is included, a large 
scale field can be generated whose structure and depen- 
dence on boundary conditions closely resemble what is 
known from mean-field aw-dynamo theory (Brandenburg 
1998). It has been suggested that nonline ar effects can 
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suppress the a-effect to insignificant values (Vainshtcin & 
jCattanco 1992| ), but as recent simulations have shown, this 
does not affect the very possibility of mean-field dynamo 
action, but merely the time-scale on which the large scale 



field is established (Brandenburg 2001). With these justi 
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fixations in mind, we adopt the mean-field approximation. 
The evolution of the magnetic field is given by 



dB 



= V x (V x B + aB - r/V x B), 



(1) 



where 77 is the turbulent magnetic diffusivity and a the 
coefficient quantifying the a-effect. We ignore anisotropies 
of a because only the toroidal diagonal component of the 
a-tensor matters in the so-called af£-regime. The antisym- 
metric part of the a-tensor corresponds to turbulent dia- 
magnetism which tends to push magnetic fields into the 
disk, and therefore just lo wers the effective win d speed. 
Furthermore, as shown by |Gabov et al. (1996 , Fig. 6), 



turbulent diamagnetism is relatively unimportant for dy- 
namo numbers less than 200, and our results mostly re- 
fer to this range. Given the substantial uncertainties in 
the mean-field transport coefficients, we decided to adopt 
the simplest possible approach. In order to minimize the 
number of parameters in our model, we also neglect the 
radial component of the velocity and thus our velocity 
field is V — (0, rCl, V z ) in cylindrical coordinates (r, ip, z) 
(see comments in Sect. 2.2 on the consequences of this 
assumption) . 

We use two approaches to solve this equation. Firstly, 
a 2D numerical simulation based on a finite-difference 
scheme is applied (Sect. 2.1). Secondly, we develop a ID 
free-decay mode expansion (Sect. 2.2). Details of these two 
complementary approaches are given below. Each of these 
two methods has its own advantages. The finite-difference 
simulation allows investigation of the nonlinear behaviour 
and applies to 2D geometry. On the other hand, the 
free-decay mode expansion gives the full eigenstructure 
and provides a tool for qualitative analysis. Furthermore, 
boundary conditions corresponding to an ideal vacuum 
can be implemented easily in the ID model, whereas our 
2D model where the disk is embedded in a halo can pro- 
vide only an approximation to a vacuum outside the disk. 
Our two approaches are therefore complementary. 

The problem is governed by three dimensionless num- 
bers which appear in the dimensionless form of Eq. ([!]). 
Let h, the half-thickness of the disk, be a characteristic 
length-scale, ceo a characteristic value of a, V z q a charac- 
teristic vertical velocity and 770 a characteristic value of 
77. We also define the rotational shear S(r) — rdfl/dr. 
The dimensionless numbers are the azimuthal magnetic 
Reynolds number 

(2) 



R 



the magnetic Reynolds number based on the a-effect 
_ a h 
Vo 

and the vertical magnetic Reynolds number 
V z0 h 



Rv — 



(3) 



(4) 



Note that ao and S depend on r and thus the magnetic 
Reynolds numbers R u and R a also depend on r. However, 



we assume later that agS is approximately constant. In 
this work, we adopt the aw-approximation for the dy- 
namo. In that case, for thin disks, the growth rates depend 
on the dynamo number T> = lZ a lZ u only and not on the 
individual values of R a and R^. However, note that the ra- 
tio of the radial to the azimuthal magnetic field |£? r |/|.E?,J 



is about -y/ 1 J?c e I / 1 -f^qj I ( Ruzmaikin et al. 1988| ). Moreover, 
Pudritz (1981 ) showed that the dynamo number T> is ap- 
proximately constant in a thin accretion disk. 

2.1. The finite-difference model 

To ensure that the magnetic field B is solenoidal, we 
evolve the induction equation in terms of the vector po- 
tential A, where B = V x A, 



dA 

— = VxB + aB- rjnoj, 



(5) 



where j = V X B/(j,q is the current density due to 
the mean magnetic field and /io is the magnetic perme- 
ability. The a-term aB together with the velocity term 
V x B drives the mean-field dynamo. We restrict our- 
selves to axisymmetric solutions of Eq. (|J). We do not 
make the thin-disk approximation in this model but we 
consider thin disks by choosing an appropriate aspect 
ratio (see below). We do not make the aw-dynamo as- 
sumption explicitly i.e. we retain the a-effect in the equa- 
tions for A r and A z , but we restrict ourselves to the au>- 
dynamo regime by our choice of magnetic Reynolds num- 
bers \Ru,\/\Ra\ = 10M0 6 . 

Our computational domain contains the disk embed- 
ded in a halo: < r < L r =2, \z\ < L z —1, where the 
disk occupies the region < r < 1.5, \z\ < 0.15, R = 1.5 
being the outer radius of the disk and h = 0.15 its half- 
thickness. This results in a disk aspect ratio oih/R = 0.1 
and the upper boundary of the domain at z ~ 7h. Using a 
resolution of 101 x 101 grid points and a uniform grid, the 
resulting mesh spacing is Sr = Sz = 0.02. Equation (||) is 
solved using a sixth order finite-difference scheme in space 
and a third order Runge-Kutta time advance scheme. 

We control the symmetry of our solutions (dipolar or 
quadrupolar) by using appropriate initial conditions. On 
the boundaries of the computational domain (i.e. far away 
from the disk) we impose the conditions that the normal 
component of A, A±, vanishes together with the normal 
derivative of the tangential components, An, 



A 



0. 



—A 

dn " 



0. 



(6) 



This is similar to the "normal-field" condition, where the 
tangential components of the magnetic field and the nor- 
mal component of the electric field E vanish, B\\ = 0, 
E±_ = 0. To verify that the boundary conditions do not 
affect our results, we have also tried perfectly conducting 
boundary conditions and found differences in the magnetic 
field structure only in the close vicinity of the boundaries. 
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The a- coefficient a(r, z) is antisymmetric about the 
midplane and vanishes outside the disk. We adopt 



a(r, z) 



ao(r) sin(7rz//i) for \z\ < h, 
for \z\ > h. 



(7) 



The radial profile of a(r, z) is smoothly cut off at r > R 
and at j"j = 0.15 where the rotational shear is very strong. 

As appropriate for accretion disks, we adopt a softened 
Keplerian angular velocity profile in r in the whole com- 
putational domain (including the halo as well as the disk), 
assuming 51 to be z-independent: 



fi(r) 



GM 



1 



-(n+l)/2n 



(8) 



where G is Newton's gravitational constant, M the mass 
of the central object, tq = 0.05 the softening radius, and 
n = 5. At r = 0, f2 vanishes as r 3 / 2 . 

The turbulent magnetic diffusivity is given by an in- 
terpolation between r^disk in the disk midplane and ?7haio 
in the halo, 



r)(r, z) = ry halo + (?7 disk - 77 halo )£(r, z). 
Here, the profile £(r, z) defines the disk, 
£,( r , z ) = g( r - r ,d)g(\z\ - z ,d), 



(9) 



(10) 



where g(x,d) is a smoothed Hcaviside step function with 
a smoothing half- width d set to 8 grid zones. So £(r, z) is 
equal to unity deep in the disk and vanishes in the halo. 

To obtain the magnetic Reynolds numbers of Eqs. @- 
(||) we choose 770 = Vdisk- 

In accretion disks, the dynamo number is approxi- 
mately constant with radius. In this model this is achieved 
by setting 



S(r)h 3 ' 



(11) 



2.2. The free-decay mode expansion 

The lowest order approximation for Eq. (Q) in a thin disk 
results from retaining derivatives with respect to z alone. 
In the aw-approximation and with a constant 77, Eq. (|TJ) 
then becomes (see Ruzmaikin et al. (1988) ) for a form with 

Vz = 0) 



dB r 
~dt 

dB v 

at 



d 

- — (aB v ) 
oz 



SB r + r] 



d 2 B, 
dz 2 



d 2 B r 
' dz 2 
d 



oz 



(y z B v ), 



(12) 
(13) 



where once again S(r) = rdtt/dr and where B r and B v are 
respectively the radial and azimuthal components of the 
magnetic field. Note that we have neglected V r , the radial 
component of the velocity in the disk, but the inclusion 



of this extra term would be straightforward (Moss et al 
2000|). 



We consider solutions with exponential time depen- 
dence 



B r (z,t) — B r (z) exp(7t), 
B v (z,t) = B v (z) exp(7i), 



(14) 
(15) 



where 7 is the complex growth rate. Measuring z in units 
of the half-thickness h of the disk, t in units of the diffusion 
time h 2 /n, the coefficient a in units of a characteristic 
value «o an d V z in units of a characteristic value V z q, the 
dimensionless form of Eq. ( |l2| ) and Eq. (13) is given by 



1 B r 



-R a ^-{aB v ) 
dz 



d 2 B r 



Rv-r{V z B r ), 
dz 



^B v — RujB r 



d 2 B, 
dz 2 



Rv-r{V z B v ), 
dz 



(16) 
(17) 



where we have kept the same notation for the dimension- 
less variables and where R a and R v are given by 
Eqs. (|), (|) andj|) with^r/o = 77. 

We expand B r and B v with respect to the orthonor- 
mal base formed by the free-decay modes (Galerkin ex- 
pansion) , 



N 

B r = a n B°' n 



N 

B<p = KB°' n 7 



(18) 



where a n and b n are complex coefficients. The free-decay 
modes are the solutions of the equation 



„,0 r>0,n u r,0.n 



(19) 



Vacuum boundary conditions are imposed at the disk sur- 
face. Wejmpose them for each free-decay mode so that 
B r and B v , and thus B r and B v , satisfy these conditions 



B' 



z=i — B' n \ z =i — 0. 



(20) 



At the midplane, the boundary conditions are connected 
to the symmetry of the mode considered. For quadrupolar 
modes we have 



dB 



0.1! 



dz 



dB / 1 



dz 



= 0, 



z=0 



and for dipolar modes we have 



rS r ' z =o — t5 U =0 — U. 



(21) 



(22) 



Thus, the normalized eigenfunctions of free decay and the 
corresponding growth rates are, for quadrupolar symme- 
try 

= %/2cos((i+n)7rz) with n = 0,1... N, (23) 
7 o = _(i +n )V, (24) 
and for dipolar symmetry 

= y/2 sin(nTrz) with n = 0,l...N, (25) 
7 ° = -n 2 n 2 . (26) 
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Note t lat the free-decay dipolar mode n = would cor- & Shukurov 1990; Poezd et al. 1993). It also allows for a 



respond to a uniform steady vertical magnetic field. The 
free-decay modes are normalized such that 



Jo 



(27) 



(except for the dipolar mode with n = 0). The 
free-decay modes are degenerate, forming doublets 
(-B°'™, 0), (0, B°' n ) decaying at the same rate 7°. 

After making this expansion, the differential equations 
become an algebraic set of equations forming an eigenvalue 
problem of 2(N + 1) equations 



N 



N 

R v J2 a n W r nm , (28) 



(7 ~ 7m) b m = Ruja m - Rv ^2 b n W% n 



(29) 



71=0 



with m = 0, 1 ... N. 

The coefficients WT m , W% m and Y nm are defined by 



W r = 



* * nm 



[ Br ,m (V z B° ,n )' dz, 
Jo 



B°> m (V z B»> n )'dz 



0,n\l , 



Y 

1 n. 



B o,m (aB o,»y dSj 



(30) 
(31) 
(32) 



where the prime denotes the derivative with respect to z. 
We normalize the eigenfunctions (Eq) such that 




N 



d Z = j: 



71=0 



= 1. (33) 



In this model, the truncation level needed depends on 
the values of R v and T>. Typically, N — 4 is enough to 
compute the eigenvalues (growth rates) of the problem 
whereas N = 9 is usually required to compute the eigen- 
functions. The precision of all the results presented here 
was checked. 

2.3. Choice of parameters 

In this section, we discuss properties that are common to 
the two models. 

The rotation of the thin disk is assumed to be close 
to Keplerian and independent of z for the finite-difference 
model (see Sect. 2.1). Only the assumption of the inde- 
pendence of f2 on z is needed in the free-decay mode ex- 
pansion. 

We use two configurations: one where the magnetic 
diffusivity is the same in the disk and in the halo around 
the disk (rjhalo = ?7disk) and one with a low conductivity in 
the halo or a vacuum around the disk. Note that assum- 
ing a vacuum or low conductivity in the halo above the 
disk is a rather good approximation for galaxies ( [Bokolofi 



connection with established theory. Note finally that the 
low-conductivity approximation (J w 0) is a good approx- 
imation for potential fields (j — 0) and hence for a par- 
ticular case of force- free fields (j x B = 0). 

The boundary conditions of the free-decay mode ex- 
pansion imply a vacuum above the disk. Technically, to 
achieve ?7haio = ?7disk in this model, we restrict the dy- 
namo action to an inner layer of the disk. The outer part 
of the disk is thus playing the role of the halo. The results 
simply need to be rescaled for a proper comparison. 

In the finite-difference model, we study the case 
??haio = ?7disk and the case of a low conductivity halo with 
?7haio/??disk = 20. Note that, with this model, to approxi- 
mate a vacuum in the halo we have to choose ?7haio/?7disk 
as high as possible. But, this choice is restricted because 
the higher this ratio, the smaller the time step. 

In this paper, the model with r^aio = ?7disk is referred 
to as the "disk+halo" model whereas the model with a 
vacuum or low conductivity above the disk is referred to 
as the "disk+ vacuum" model. 

In both models, the profile for a is taken to be 



ao sin 7r 



for 



< h 



(34) 



(see Eq. (Q) for details concerning the finite-difference 
model), except in Fig. [j] where we take a linear depen- 
dence of a with z to compare with previous work. The 
profile for V z is linear in z, 



V z = V zQ ~. 
h 



(35) 



The accuracy of the models depends on the number of 
grid points in r and z for the finite-difference model and on 
the level of truncation for the free-decay mode expansion. 

In this work, the two control parameters are the dy- 
namo number T> and the vertical magnetic Reynolds num- 
ber R v . For reference, we give here possible values for these 
two parameters. While dynamo numbers are rather low in 
galactic disks, \T>\ < 10 (Ruzmaikin et al. 1988), in ac- 
cretion disks, they can be as high as \D\ ~ 10 2 — 10 4 , 
depending on the Mach number of the turbulence, see e.g. 
Pudritz (198l[ ) . The value of R v obviously depends on the 
origin of the vertical outflow. With magnetic buoyancy, 
a rough estimate can be obtained assuming that the ver- 
tical motions occur at the Alfven speed. Estimating the 
Alfven speed from magnetic equipartition, one gets a ver- 
tical magnetic Reynolds number R v of order unity. 

In the following, the figures obtained with the finite- 
difference model are labelled with "FD" and the ones 
obtained with the free-decay mode (Galerkin) expansion 
with "GE". 



3. Linear results for R r 







In this section, we present results with no vertical velocity 
in the disk (R v — 0). Such a study is an interesting test for 
the quality of the results of the free-decay mode expansion. 
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It also gives a point of comparison between the numerical 
simulation and the free-decay mode expansion. 

We first compare the results of the free-decay mode 
expansion (with R v — 0) with the results published by 
Stepinski fc Levy (199l| ) . We thus take a linear dependence 
of a on z, a — aoz/h. The calculated growth rates are 
plotted in Fig. |[ The agreement is good but a discrepancy 
is observed at large \T>\ where the increase of the growth 
rate with |X>j is faster according to Stepinski & Levy than 
in our free-decay mode expansion for both quadrupolar 
and dipolar symmetry (T> < 0). This discrepancy seems to 
be due to an insufficient number of grid points in Stepinski 
& Levy's model. Indeed, if we increase the truncation level 
of the free-decay mode expansion, no change is observed 
in the results whereas if we decrease the truncation level, 
we recover the results of Stepinski & Levy. Note that the 
results discussed in this paper are for \T>\ < 400 where the 
agreement is excellent. 

We now want to compare the finite-difference simula- 
tion with the free-decay mode expansion. In particular, we 
want to test our method to mimic, in the case of the free- 
decay mode expansion, a model for ?/haio = f?disk starting 
with the "disk+vacuum" model (see Sect. 2). For an in- 
finitely extended halo with ryhaio = ?7disk, the growth rates 
of dipolar and quadrupolar modes should be exchanged 
when T> reverses sign, i.e. the graph of -Re (7) as a func- 
tion of T> should be symmetric with respect to the vertical 
axis ( Proctor 1977 ). The finite-difference simulation shows 
this symmetry with a good precision (Fig. |^a). The results 
of the free-decay mode expansion (Fig. |2Jb) are less sym- 
metric. This is not surprising since the ratio between the 
height of the computational domain and the height of the 
disk is about 7 in the finite-difference simulation whereas 
it is only 2 in the free-decay mode expansion. 

We now discuss our results for the "disk+halo" and 
"disk+vacuum" configuration. Let us first consider the 
"disk+halo" configuration (Figs. ||a and b). For moderate 
dynamo number |X>| < 300, the leading mode is quadrupo- 
lar and non-oscillatory if T> < (R a > 0) and it is dipolar 
and non-oscillatory if T> > (R a < 0). In the following, 
we thus study the effect of vertical motions for V < in 
quadrupolar symmetry and for T> > in dipolar symme- 
try. Above \T>\ ~ 300, the symmetry of the leading mode 
changes and it becomes oscillatory. The critical dynamo 
number is |2? C rit| ~ 5. 

Let us now consider the "disk+vacuum" configuration 
(Figs. |2|c and d). The results obtained with the two ap- 
proaches are rather different. With the free-decay mode 
expansion, the dipolar mode is always oscillatory and sub- 
critical for \T>\ < 400. On the other hand, with the finite- 
difference simulation, the dipolar mode becomes supercrit- 
ical for V > 50 and is non-oscillatory for all V > 0. This is 
because the dominant dipolar mode in the finite-difference 
simulation is the "forgotten mode" i.e. a dipolar mode cor- 
responding to rt = in Eqs. ( p5|) and (26). The "forgotten 
mode" cannot be obtained with our ID calculation (free- 
decay mode expansion). However, we reproduced the cor- 
responding profile 7(2?) with a modified ID model which 




100 200 300 400 500 600 700 800 900 1000 

\D\ 

Fig. 1: Real part of the growth rate of the magnetic field 
as a function of \T>\, with a — aoz/h, obtained with the 
free-decay-mode expansion with N > 9. Asterisks denote 
non-oscillatory solutions whereas diamonds denote oscilla- 
tory solutions. Solid lines are for quadrupolar modes and 
dotted lines for dipolar modes. Signs above the curves in- 
dicate the sign of D. 



takes into account the finite aspect ratio of the disk. The 
model is based on the assumption of separability of r- and 
z-dependence and yields equations similar to the thin-disk 
equations (|lj) , ( [l3| ) , but with additional correction terms 
of order (h/r) 2 . Thus, the dipolar modes of Fig. ||c and 
Fig. ||d are of a different nature. As for the quadrupolar 
mode, the results obtained with the two approaches are 
not so different. The branch for T> < is rather similar 
in both cases with a critical dynamo number > —50. The 
branch for T> > is all subcritical with the finite-difference 
simulation whereas it becomes supercritical for T> ~ 170 
with the free-decay mode expansion. However, we have 
checked that, in the finite-difference model, if we take a 
steeper transition between ?7disk and ryhaio, the branch for 
T> > moves upwards, resembling the one for the free- 
decay mode expansion. 

The qualitative behaviour for the case of low conduc- 
tivity in the halo is in many respects similar to that of 
the "disk+halo" configuration (Figs. ||a and c). But there 
are quantitative differences. There is no intersection point 
where the symmetry of the leading mode would change 
for \T>\ < 700. The diagram does not look symmetric any- 
more. 

The results of the finite-difference model presented 
here are in agreement with the accretion disk models of 
Rudigcr ct al. (1995| ) and |von Rekowski et al. (20"oo| ). Note, 
however, that these models have an inner radius and that 
at this radius a perfectly conducting boundary condition 
is used, which leads to a vertical current layer on the 




Fig. 2: a-d Real part of the growth rate of the magnetic field as a function of —T> for R v — 0. Asterisks denote non- 
oscillatory solutions and diamonds denote oscillatory solutions. Solid lines indicate quadrupolar modes and dashed 
lines indicate dipolar modes, a and c are for the finite-difference model (FD) and b and d are for the free-decay mode 
expansion model (GE) with a and b for the "disk+halo" model and c and d for the "disk+vacuum" model. Note the 
approximate symmetry of graphs a and b with respect to the vertical axis. Note also that the dipolar mode in c is 
the "forgotten mode" , not seen in d. 



boundary and potentially lowers the excitation threshold 
for magnetic fields. 



4. Linear behaviour with vertical velocities 

In this section, we study the linear behaviour of the aui- 
dynamo in the presence of vertical velocity (R v ^ 0). As 
discussed in Sect. 3, we present results for quadrupolar 
symmetry with T> < (Sect. 4.1) and results for dipolar 
symmetry with V > (Sect. 4.2.). An explanation of these 
results is given in Sect. 4.3. 



4.1. Quadrupolar modes 

We consider quadrupolar modes with T> < 0. For suffi- 
ciently low dynamo numbers (|X>| < 50), the growth rate 
decreases with R v whereas at larger dynamo numbers, the 
growth rate increases for small R v and then decreases at 
large R v , forming a maximum (Fig. [3]). This maximum is a 
robust feature, the existence of which does not seem to de- 
pend on the boundary conditions or on the model. Indeed, 
we find this maximum with the finite-difference simulation 
(Figs. 3a and c) and with the free-decay mode expansion 
(Figs. 3b and d) in both configurations, "disk+halo" and 
"disk+vacuum" (see Sect. 2.3 for the definition of these 
configurations). The boundary conditions at the surfaces 
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of the disk affect the position of the maximum in R v and 
the dynamo number above which it occurs. In particu- 
lar, the smaller the magnetic diffusivity of the halo, the 
smaller this dynamo number. Nevertheless, the required 
values for T> to get a maximum are larger than the typical 
values in galaxies. When T> is increased, the maximum is 
more pronounced and its position is shifted towards larger 
Rv. Note that the values of R v at which the maximum oc- 
curs are between 1 and 10. They thus fit within the range 
of values estimated for magnetic buoyancy. Finally, in the 
range of dynamo numbers studied here, \T>\ < 300, the 
vertical velocity does not change the symmetry or tempo- 
ral behaviour of the dominant mode. 

We now discuss the magnetic field structure. Figure |] 
shows the eigenfunctions for B r and B v for T> = —200 and 
for different R v in the "disk+halo" configuration. With 
this choice of the dynamo number, a maximum in the 
curve j(Rv) occurs around R v = 7.5. An obvious feature 
is that the vertical velocity pushes the zeros of B r and B v 
outwards. For the dynamo to work, it is essential that B r 
changes sign in the disk since only then can magnetic flux 
of the opposite sign to that in the midplane leave through 
the disk surface ( Ruzmaikin et al. 198§| ). When R v takes 
large values, the radial component of the magnetic field B r 
has no zero inside the disk and thus the dynamo cannot 
work (7 < 0). This is the technical reason for the decrease 
of 7 with R v . The physical reason for it is investigated 
in Sect. 4.3. The increase of the growth rate at small R v 
will be explained in Sect. 4.3. We note however that when 
we decrease the truncation level of the free-decay mode 
expansion down to N — (one doublet), the maximum 
disappears. 




Fig. 6: Real part of the growth rate of the magnetic field 
as a function of R v for the dominant dipolar mode and 
for different dynamo numbers in the "disk+halo" model. 
Asterisks denote non-oscillatory modes whereas diamonds 
denote oscillatory ones. 




The magnetic field configuration in 2D geometry is 
shown in Fig. ||. It shows both the effect of increasing 
the dynamo number and the effect of increasing the ver- 
tical magnetic Reynolds number. When one decreases the 
dynamo number from T> = —50 to I? = —200, in the ab- 
sence of any vertical velocity, the magnetic field becomes 
more concentrated at the outer edge of the disk (Figs. ||a 
and |5|b) . If the vertical magnetic Reynolds number is then 
increased from R v — 0, the magnetic field is radially redis- 
tributed over the whole disk (Fig. |B|c) and then advected 
towards the halo (Fig. ||d). From panel (b) to (d), the 
growth rate increases as shown by the curve for T> = —200 
in Fig. Ha. Thus, the magnetic field diffusion becomes less 
efficient as the vertical scale of the magnetic field increases 
with R v . This leads to an increase of the growth rate up 
towards its maximum at R v — 7.5 (Fig. gd). Note that 
in the halo the poloidal magnetic field lines become more 
and more aligned with lines of constant angular velocity as 
R v increases. Note also that the mode structure of Fig. ||d 
(where Re(7) is maximum for T> = —200) is very similar 
to the mode structure for — T> = 50, R v = 3, where that 
mode has attained its maximum growth rate. 



Fig. 7: a and b. Magnetic field configuration obtained us- 
ing the finite-difference model with the linear "disk+halo" 
configuration and for the dominant non-oscillatory dipolar 
mode. The parameters are a T> = 50, R v = 0, b T> = 50, 
R v — 0.75. The notation is the same as in Fig. B. 

4.2. Dipolar modes 

We now study the linear behaviour with vertical velocity 
for positive T> where the dipolar non-oscillatory mode is 
dominant. Figure ^ shows that there is no maximum in the 
growth rates; at R v of order 1 or even less the dynamo is 
switched off. 

The magnetic field structure in 2D geometry is shown 
in Fig ^. The magnetic field structures for V = 50 and 
R v — and 0.75 show that the field is advected vertically 
outwards and the poloidal field lines are aligned with the 
fl contours lines very quickly. 

4.3. The origin of the dynamo enhancement 

In this section, we explain why a vertical velocity can in- 
crease the growth rate of the dynamo. Therefore, we plot 
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Fig. 3: a-d Real part of the growth rate of the magnetic field as a function of R v for the dominant quadrupolar mode 
and for different dynamo numbers. All modes are non-oscillatory, a and c are for the finite-difference model (FD) and 
b and d are for the free-decay mode expansion model (GE) with a and b for the "disk+halo" model and c and d for 
the "disk+ vacuum" model. The models of panels a to d correspond to the ones in Fig. 0. 



the different terms of Eq. ([l6]) and Eq. ( |r7|) obtained with 
the free-decay mode expansion. To avoid any ambiguity 
with the choice of normalization, in particular when com- 
paring results at different R v , we rewrite the equations in 
the following way. We multiply Eq. ( |l6| ) and Eq. ([l7]) by 
B r and B v respectively, and we integrate over the disk 
thickness. We thus obtain 



7 = la + 7diff + 7stret> 



7 = 7^ 



with 



■ Tdiff " 



/strct J 



^B r {aB v )'dz 
Jo B*6z 



(36) 
(37) 



(38) 



7 di ff 



/strct 



j B r B^dz 
Jo B^B'^dz = _JoiK 

Jo B r „ (V z B r>v )'dz r v /o 1 v;_B 2 r v dz 

' 2 J^ dz 



r,<p ) 



-Rv 



Jo B r,<p dz 



(39) 
(40) 
(41) 



where we have used a vacuum boundary condition at z = h 
to rewrite 7strct, 7diff- This form of the equations allows 
us to investigate the different physical contributions to 
the growth rate 7. Note that Eqs. (|3^) and ( |37| ) are not 
explicitly time-dependent. The terms 7^ and 7^ are the 
contributions to the growth rate 7 from the a-effect and 
w-effect. The terms 7^ are the contributions from dissi- 
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Fig. 4: a-d Vertical structure of the non-oscillatory quadrupolar magnetic field for T> = —200. Panels a and c are for 
the finite-difference model (FD) while panels b and d are for the free-decay mode expansion (GE). Panels a and b 
give the radial structure and c and d give the azimuthal structure of the magnetic field. In the finite-difference model, 
eigenfunctions are taken at r = 1 and are normalized such that max \ B V \ = 1 for each R v . Note that in the free-decay 
mode expansion, we plot B r and B v , i.e. the time-independent part of the magnetic field. 



pation which are always negative. Finally, the terms 7j^t 
are the contributions to 7 from stretching the magnetic 
field. The profile for V z is always antisymmetric with re- 
spect to the disk midplane and Vg > for z > 0. Thus, 
these stretching terms 7^ret are always negative. In par- 
ticular, with a linear profile for V z , they are both equal to 
-Rv/2. 

The general picture is as follows. The positive contri- 
bution to 7 comes from 7^ and 7^ whereas dissipation has 
a negative contribution. In the presence of a vertical ve- 
locity, an additional negative contribution comes from the 
stretching of the magnetic field %if ct - But, at the same 
time, the stretching reduces dissipation by increasing the 
scale of the magnetic field (see Eq. (fib])). Thus, depend- 



ing on the relative values of the reduction of dissipation 
and of the negative stretching term, the growth rate is 
enhanced or not (a maximum occurs or not). Obviously, 
when the magnetic field is completely stretched, dissipa- 
tion cannot be reduced anymore. However, this does not 
necessarily correspond to the maximum of the growth rate 
as the terms 7^ and 7^ evolve themselves with R v (which 
quantifies the vertical velocity). These different quantities 
are plotted in Fig. | in a case where a maximum occurs 
and in a case with no maximum. 

We note that dissipation increases again after being re- 
duced by the stretching of the magnetic field in the vertical 
direction. This increased dissipation is due to the forma- 
tion of boundary layers near z = h. We also note that the 
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Fig. 5: a-d Magnetic field configuration obtained using the finite-difference model with the linear "disk+halo" config- 
uration and for the dominant non-oscillatory quadrupolar mode. Contours represent poloidal magnetic field lines with 
solid lines being clockwise and dotted lines being counter-clockwise. The greyscales are the intensity and direction 
of the azimuthal magnetic field; light grey (dark grey) signifies B v > (B v < 0). The parameters are a T> = —50, 
R v = 0, b V = -200, R v = 0, c V = -200, R v = 0.75, dV = -200, R v = 7.5 (where Re(j) is maximum). Comparing 
a and b shows the effect of increasing the dynamo number, b-d show the effect of increasing R v . 



term 7^ decreases whereas 7^ increases. The evolution of 
these terms is linked to the disappearance of the positive 
area of B r (with the convention of sign of Fig. ||). Finally, 
note that at large R v , the growth rate 7 always decreases 
as the stretching terms 7^ t become dominant. 

In short, the growth rate can increase with R v , thus 
forming a maximum, because dissipation is reduced by the 
stretching of the magnetic field. 



5. Nonlinear behaviour with vertical velocities 

The nonlinear behaviour has only been investigated with 
the finite-difference model. We restrict ourselves here to 
the model "disk+halo" (?yh a io = ?7disk) and to negative 
dynamo numbers considering only quadrupolar symmetry 
which, for R v = 0, is the dominant symmetry for moder- 
ate values of V < 0. Solutions are all saturated and non- 
oscillatory, as in the corresponding linear modes. We study 
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radial equation, £> = -200 radial equation, D=-10 




Fig. 8: a-d Contributions of the different terms in Eqs. ([36]) and ( |37| ) to the growth rate 7 in the "disk+ vacuum" 
model. Lines are labelled with the subscripts used for the different terms in Eqs. (^) and |37]). a and c represent the 
radial and azimuthal terms for T> — —200, for which a maximum occurs, b and d represent the radial and azimuthal 
terms for D = —10, for which no maximum occurs. 



three different cases of nonlinearity: magnetic buoyancy, 
a-quenching and a combination of both. For a compari- 
son of nonlinear solutions with a different parametrization 
of magnetic buoyancy see earlier works by Torkelsson & 
Brand enburg (1994a| ,|bl) . 



5.1. Magnetic buoyancy 

We parameterize the effect of magnetic buoyancy by as- 
suming that the vertical velocity depends on the horizon- 
tal magnetic field, 



V z (z) = V z0 -jB?+Bl/B c , 



where B cq is the equipartition magnetic field £> oq = 
y/p/IoCsMaturb, P is the density c s the sound speed and 
Maturb a turbulent Mach number, treated here as a free 
parameter. The profiles of density p and sound speed c s 
are taken from a simple model of a cool accretion disk 
embedded in a hot corona (Brandenburg et al. 2000; see 
also [Doblcr et al. 1999|). 



We define an effective vertical magnetic Reynolds num- 



ber 



Ri 



oil' ) 



R v JB* + Bl/B 



disk 



(43) 



where (-)disk denotes averaging over the disk. Fig. ^ shows 
the magnetic energy and i?i eff ' as a function of time for 
(42) d = —200 and R v = 7.5. This value of R v corresponds to 
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Fig. 9: Magnetic energy and effective vertical magnetic 
Reynolds number R v cS ^ as a function of time for the 
"disk+halo" model with magnetic buoyancy (cf. Eq. fl42])). 
T> = —200 and R v — 7.5. Super-exponential growth oc- 
curs before saturation sets in. The mode is quadrupolar 
and non-oscillatory. The initial magnetic field strength is 
-Binit = 10 -5 . The steep increase for t < 2 is due to the 
generation of azimuthal field from the initial large scale 
poloidal one at the chosen high R u . 

the location of the maximum of the growth rate j(Rv) in 
the linear model (cf. Fig. upper curve). Starting from 
an initial magnetic field strength B m n = 10~ J , B grows 
and therefore R v eS ^ increases; so the growth rate roughly 
follows the upper curve in Fig. ||a. Thus, as long as Rv 
is approximately less than the position of the maximum 
of 7, i.e. the horizontal magnetic field is smaller than the 
equipartition field, 7 increases with time, which results 
in the super-exponential growth of the magnetic energy 
seen for 10 < t < 28 (Fig. |). After t « 28, when the 
horizontal field is comparable with the equipartition field, 
7 decreases and eventually the magnetic field approaches 
its saturation level. This level is reached when R v eS ^ ~ 
R v * « 21.5, where R v * is the zero of ~f{R v ) in Fig. ||a, 
upper curve. 

In this subsection we assume that the only nonlinear- 
ity is due to magnetic buoyancy and thus appears in the 
vertical velocity term only. In particular, we do not con- 
sider here any dependence of a or 77 on B. In this case 
the dimensionless form of the induction equation in the 
steady state reads 



£B + Vx 



R, 



Bl+Bl/B C A uxB 



= 0, 



(44) 



where C is a linear operator and u is a vector function of 
position, u — (0, 0, z/h) in the case considered here. Let b 
be the solution of Eq. ( p4| ) for R v — 1 . Then, for arbitrary 
R v , the solution is b/R v . From this, it follows that the 
saturation value of magnetic energy is 



Et 



_ i^s at 

Rl mag 



1 



=1* RV 



(45) 




Fig. 10: Magnetic energy as a function of R v for the 
"disk+halo" model with a-quenching. All modes are 
quadrupolar and non-oscillatory. 

and the structure of the saturated magnetic field is in- 
dependent of R v . This is confirmed with our numerical 
model. 

The magnetic field is extended into the whole halo and 
almost completely stretched. Therefore, adding a wind 
term such that 



v.(*) = (vjry 



Bl 



Bl/B 



eq 



v.. 



wind 
zO 



(46) 



and fixing R^ = R v cS) while increasing Rf nd = R v , 
cannot lead to enhancement of magnetic energy. The wind 
is carrying magnetic field away and -Emag is decreasing. 

5.2. a-quenching 

We now parameterize the back-reaction of the magnetic 
field on the turbulent a-effect by making a depend on B, 



a(r, z) = cv (r) sin y~ 



■B*/Bl 



(47) 



where B cq is again the equipartition magnetic field. 



Figure 10 shows the magnetic energy as a function of R v 
for the three dynamo numbers T> = —50, T> = —100 and 
T> = —200. From this we can see the following features: 



- EL 



has a maximum as a function of R v in the non- 



linear regime; 

— but, this maximum is shifted towards much smaller R v 
(of order 0.1), compared to that of the growth rate in 
the linear regime, and it is also less pronounced; 

— the value of R v where the maximum in -E mag occurs, 
is almost independent of T>\ 

— the magnetic energy scales roughly with the dynamo 
number, i^mag oc T>; 

— all modes remain non-oscillatory when increasing R v . 



z 0.0- 
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Fig. 11: a-c Magnetic field configuration for the "disk+halo" model with a-quenching: the effect of vertical velocity. 
Poloidal field (lines) and azimuthal field (grey scales). aT> — —50, R v = 0; b T> = —50, R v = 3; c T> = —50, R v = 7.5. 
a-c following the lower curve in Fig. ffOl All modes are quadrupolar and non-oscillatory. 



The effect of vertical velocity on the structure of the mag- 
netic field for the model with a-quenching can be seen in 
Fig. [ll] where it is shown for V = —50. But looking at the 
magnetic field configuration for T> = —100 and T> = —200, 
one finds that the field structure changes only weakly with 
the dynamo number. 

5.3. a-quenching together with magnetic buoyancy 



In this section, we combine a-quenching, Eq. (|47|), with 
magnetic buoyancy, Eq. ([42|), and compare the results with 
the case of a-quenching with wind. In Fig. 12, magnetic 
energy is plotted against R v cS ^ for these two models and 
for V = —50. Evidently, the two curves behave similarly, 
supporting the notion that wind and magnetic buoyancy 
have similar effects on the dynamo. For R^ S ^ <C 1, the 
field structure is solely determined by a-quenching and 



thus -E mag is independent of R\, 



(off) 



For R v cS) 



» 1, 



the other hand, a-quenchmg becomes unimportant and 
-Bmag cc (i?i cff ')~ 2 , as discussed in Sect. 5.1. The largest 
disagK ement between the two curves occurs for intermedi- 
ate val res of r[ cS \ i.e. where both a-quenching and mag- 
netic buoyancy are important, and where the field has not 
yet reached its asymptotic shape (Fig. [Tl]c) . 

Perfect agreement between the two curves cannot be 
expected, because the averaging in Eq. (^3|) was somewhat 
arbitrary; the agreement could perhaps be improved by 
choosing a different definition of Rv . The fact that with 
magnetic buoyancy the magnetic energies are larger for a 
given value of R v cti ^ can therefore not be interpreted as a 
relative enhancement of dynamo action. 



6. Conclusions 

Disk dynamos have been studied intensively over the past 
20 years, but only recently the effects of outflows and 
winds have been included (e.g. Campbell 1999; Doblcr 
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Fig. 12: [Magnetic energy as a function of Rv defined in 
Eq. ( f43| for the models with a-quenching with wind and 
a-quenctiing combined with magnetic buoyancy. f?haio = 
f^disk an I T> = —50. All modes are quadrupolar and non- 
oscillatc ry. With a-quenching alone, i?i° ff ^ = R v . 



et al. 1999; |Brandenburg et al. 200C| ). It was immediately 
clear that these changes in the boundary conditions af- 
fect the structure of the solutions. The present work is 
an attempt to provide a more thorough survey of solu- 
tions with wind for different values and sign of the dy- 
namo number. One of the effects that we have highlighted 
is the enhancement of the kinematic dynamo growth rate 
for intermediate strengths of the wind. In the nonlinear 
case this corresponds to a phase shortly before saturation 
where the growth becomes super-exponential. 

In the present work we have deliberately isolated the 
dynamo problem from the wind problem. In reality the 
two aspects together constitute a more complicated self- 
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consistent dynamo problem that needs to be looked at in 
mo re detail in future. One ste p in that direction was taken 
by Brandenburg et al. (2000 ) where the feedback on the 
flow from the Lorentz force was taken fully into account. 
A more formidable task would be to obtain self-consistent 
dynamo action in three dimensions without assuming an 
a-effect, thus obtaining dynamo action and outflows in a 
self-consistent manner from first principles. In the context 
of accretion disks, this would require simulating turbu- 
lent dynamo action, where the turbulence itself is driven 



by dynamo-generated magnetic fields ( Brandenburg et al 
1995|)~|This would be of great interest not only from the 
astrophysical point of view, but also from that of dynamo 
theory in general, because, as was recently realized, the 
catastrophic quenching of the standard aw-dynamo could 
be alleviated in the presence of open boundary conditions 
which may substantially improve the working of the dy- 



namo (Blackman & Field 2000). This is simply because 
the a-effect generates helical fields, but at the same time 
magnetic helicity is conserved in a closed domain. Thus, 
being able to get rid of helical fields through open bound- 



aries may be vital for the dynamo (Brandenburg & Dobler 
200l|)7}<urther more, the production of fields at large scales 
could be facilitated by the production - and subsequent 
loss through boundaries - of small scale fields with oppo- 
site sign of magnetic helicity. This latter aspect could not 
be addressed with mean-field models and would require a 
direct simulation of accretion disk turbulence in a global 
model. 
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